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Abstract. - The gravitational hydrodynamics of the primordial plasma with neutrino hot dark 
matter is considered as a challenge to the bottom-up cold dark matter paradigm. Viscosity and 
turbulence induce a top-down fragmentation scenario before and at decoupling. The first step is 
the creation of voids in the plasma, which expand to 37 Mpc on the average now. The remaining 
matter clumps turn into galaxy clusters. Turbulence produced at expanding void boundaries 
causes a linear morphology of 3 kpc fragmenting protogalaxies along vortex lines. At decoupling 
galaxies and proto-globular star clusters arise; the latter constitute the galactic dark matter halos 
and consist themselves of earth-mass H-He planets. Frozen planets are observed in microlensing 
and white-dwarf-heated ones in planetary nebulae. The approach explains the Tully-Fisher and 
Faber- Jackson relations, and cosmic microwave temperature fluctuations of micro-Kelvins. 
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Introduction. — Structure formation in the Universe 
starts in the plasma of protons, electrons, He atoms and 
neutrinos, that exists up to some 400,000 yr after the Big 
Bang, the time of decoupling (dc) of photons from matter 
(last scattering (L) or recombination). Then the plasma 
transforms to a neutral gas of H and 24% by weight 4 He, 
with the neutrinos remaining free streaming. As this oc- 
curs at about four thousand degrees Kelvin, a moderate 
plasma temperature, we shall seek an explanation in terms 
of plasma physics and gravitational hydrodynamics alone. 
This embodies a return to the top-down scenario of large 
scale structure formation. 

Currently it is assumed that cold dark matter (CDM) 
also exists and, clustered before decoupling, has set seeds 
for baryon condensation. The so-called the concordance 
or ACDM model involves also a cosmological constant or 
dark energy. It describes a hierarchical bottom-up ap- 
proach to structure formation, stars first, then galaxies, 
clusters, and, finally, voids. 

But observations of dense clumps of ancient small stars 
in old globular clusters (OGCs) in all galaxies contradict 
the ACDM predictions that star formation should begin 



only after about 300 million years of dark ages and that 
the first stars should be 100-1000 Mq population III su- 
perstars. OGCs do not spin rapidly so they cannot be 
condensations, and their small stars imply gentle flows in- 
consistent with superstars. Other difficulties are posed by 
empty supervoids with size up to 300 Mpc reported from 
radio telescope measurements [1], dwarf galaxies with a lot 
of dark matter [2] and a preferred axis of evil spin direc- 
tion (AE) that appears at scales extending to 1.5 Gpc, 
a tenth of the horizon scale [3]. Nearly every month 
new observations arise that pose further challenges to the 
ACDM paradigm: Correlations in galaxy structures [4]; 
absence of baryon acoustic oscillations in galaxy-galaxy 
correlations [5] ; galaxies formed already when the universe 
was 4-5 billion years old [6] ; dwarf satellites that swarm 
our own galaxy just like its stars [7]. 

The recent conclusion by one of us that dark matter 
particles must have mass of a few eV and probably are 1.5 
eV neutrinos [8], means that dark matter is hot (HDM), 
urging once more for an explanation of structure formation 
from baryons alone, without a cold dark matter trigger. 

We shall discuss such baryonic clustering due to a vis- 
cous instability in the plasma, overlooked by the currently 
popular linear models of structure formation. CDM is 
assumed not to exist, while HDM, though initially impor- 
tant to maintain the homogeneity of the plasma, has no 
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role in the structure formation. Central in our discussion 
will be the huge plasma viscosity v ~ 5 10 27 m 2 s _1 arising 
from photons that scatter from free electrons. This makes 
the plasma increasingly viscous, while it is also expanding 
with space. At some age before the decoupling an insta- 
bility creates the first structures, proto-voids and proto- 
galaxy-clusters. At decoupling the viscosity drops to hot 
gas values ~ 10 13 m 2 s _1 , which creates further structures 
at the Jeans scale and at the new, small viscous scale. 

The plan of this Letter is to review modern theory of 
gravitational hydrodynamical structure formation, to eval- 
uate the estimates for various fragmentation scales within 
Friedman cosmology, and to compare with observations. 

Hydrodynamics. — The description of gravitational 
structure formation starts with Jeans 1902. He proposes 
that scales for gravitational condensation of a uniform 
fluid of density p must be larger than the Jeans acous- 
tic scale Lj = Vs/(pG) 1 ^ 2 , where Vs ~ c/V3 is the sound 
speed of the plasma and G = 6.67 10 _11 m 3 /kgs 2 is New- 
ton's constant. As the plasma Jeans scale is always larger 
than the horizon scale of causal connection, this forbids 
gravitational structure formation before decoupling. The 
Jeans criterion reflects a linear gravitational instability 
from acoustics, but it neglects the fact that self gravi- 
tational instability of a gas is absolute [9]. All density 
variations will grow or decrease unless prevented by the 
viscous forces, turbulent forces and diffusion effects. 

Conservation of the specific momentum in a fluid is ex- 
pressed by the Navier-Stokes equation, 

dv -> — 

— =X7B + VXlJ + F mscous + Father, (1) 

averaged over system control volumes exceeding the mo- 
mentum collision length scale. B is the Bernoulli group 
of mechanical energy terms B = p/p + \v 2 + Iw and the 
viscous force is F viscous = v s V 2 v + {\v s + v b )V ■ (V • v), 
with kinematic shear viscosity v s = rj/p and bulk viscosity 
v b = CI P, while other fluid forces may arise. The inertial- 
vortex force per unit mass vxlj, with uj = Vx v, produces 
turbulence if it dominates the other forces; for example, 
R e = | if x ti\ / \F v i S cous\ is the Reynolds number. A large 
viscosity corresponds to a small Reynolds number, with 
universal critical value Rg ~ 25 — 100. For adiabatic flows 
the "lost work" term Iw due to frictional losses is negligible 
so the enthalpy p/p decreases or increases to compensate 
for changes in the kinetic energy per unit mass \v 2 . 

The turbulence problem is put in a new perspective by 
Gibson [10]. Universal similarity laws are explained in 
terms of the inertial vortex forces vx uj. From Eq. (1), 
turbulence is defined as an eddy-like state of fluid motion 
where the inertial- vortex forces of the flow are larger than 
any other forces that tend to damp the eddies out [11,12]. 
By this definition, irrotational flows are non-turbulent. All 
turbulence then cascades from small scales to large be- 
cause vorticity is produced at small scales and adjacent 
eddies with the same spin induce inertial vortex forces 



that cause the eddies to merge and form bigger structures. 
Thus, turbulent motions and energy always cascade from 
small scales to large, contrary to standard turbulence the- 
ories that include in turbulence also irrotational flows and 
motions dominated by other forces. 

Gravitohydrodynamics (GHD). — In his approach 
with hydrodynamic and diffusive modelling, Gibson 1996 
[10] derives several gravitational Schwarz length scales for 
structure formation by Kolmogorian dimensional analy- 
sis. With T g = 1/^fpG the gravitational free fall time 
in the Jeans length Lj = Vsr g , there occur first the vis- 
cous length Lsv — T g\/iv, wnere v the kinematic viscos- 
ity and 7 the rate of the strain, i. e., the magnitude of 
eij = ^(dvi/dxj + dvj/dxi). Second, there is the turbu- 
lent length Lst = (er 3 ) 1 / 2 , where e is the rate of energy 
dissipation per unit mass, and third, the diffusive length 
Lsd = yjDr g , where D is the diffusion coefficient. Within 
the acoustic horizon scale, dfj ~ ct structures can form at 
scale L if djj > L > max(Xsv, Lst, Lsd)- 

We shall evaluate these scales within the flat Friedman 
metric ds 2 = c 2 dt 2 — a 2 (t)dr 2 . The Friedman equation for 
baryonic and neutrino matter reads 

^ sn <°> = nA + ^- + (2) 

where fi„ becomes a-dependent beyond the Compton tem- 
perature ~ 17,000 K. With dt = da/H ay/n(a) the age 
is J Q dt, while the 'angular' distance from us to an object 

at redshift z = 1/a— 1 reads <1a{z) = \cj (1 + z)\ C° dt/a. 
We adopt the Hubble parameter h = i?o/[100km/sMpc] 
=0.744 favored in Ref. [8], so that m v = 2 z / A {G F ) 1 l 2 m 2 e = 
1.4998 eV and Q u = 0.111//r 3 / 2 = 0.173. For baryons we 
take fl B = 0.02265//J 2 = 0.0409 from WMAP5 [13], while 
for photons fi 7 = 2.47 • lO" 5 /^ 2 = 4 - 4 6 ■ 10~ 5 . Finally, 
Ha = 1 — Slv — &b — ^7 = 0.786 assures a flat space. 

Viscous instability in the plasma. — GHD starts 
with acknowledging the importance of the photon viscos- 
ity. Because it strongly increases in time, already before 
decoupling the plasma becomes too viscous to follow the 
expansion of space [10]. Thus a gravitational instability 
occurs, that tears the plasma apart at density minima, 
thus creating voids. Cosmic (super) voids surround us at 
any distance and the furthest observable ones are located 
at the decoupling redshift. Presently, voids have a 20 times 
undcr-density with respect to the critical density. In be- 
tween voids the galaxy clusters are located on "pancakes" 
that join in superclusters. 

The shear viscosity v = v s reads for fc^T <C m e c 2 [14] 

v = V_ = J_ 5m 2 C(3)(fc B T) 4 
Pb Pb 9TT 3 h 5 c 2 a 2 m n e ' 

with m e the electron mass, a cm = 1/137 the fine structure 
constant and n e = 0.76pb/itin the electron density. With 
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n e <~ ps ~ T 3 , v increases as 1/T 2 . At WMAP5 values for 
decoupling it reaches the huge value 5.85 10 27 m 2 /s, while 
the bulk viscosity ~ 10 14 m 2 /s is much smaller. In the 
plasma the acoustic speed Vg — c/ \J3{\ + R z ) with R z = 
3f2s/4f2 7 (l + z) [15] sets the acoustic horizon scale 



d a H c (z) 



t(z) 
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(4) 



At WMAP5 decoupling it takes the value d^ = 128 kpc 
but, estimating 7 = Vs/d^, the viscous length L sv = 
{vVs/GpBd ^) 1 / 2 is then only 76 kpc, showing that an 
instability has occurred. This causes an often over- 
looked baryonic structure formation in the plasma [10]. 
The crossover of Lsv d^ occurs when d^ — 

(V s v/Gp B ) 1/3 - This happens at z vf = 5120, where d<fi = 
7.3 kpc is the initial void scale. It expands by a factor 
1 + z v { to become 37 Mpc now, a typical void size, smaller 
than the supervoids of 50-300 Mpc observed by radio tele- 
scopes. ACDM models predict such voids formed last and 
full of debris, rather than first and empty as observed [1]. 
Foreground voids will play a role in the cosmic microwave 
background (CMB) structure at large angles [16], espe- 
cially due to their neutrino depiction at z <~ 28 [8] . Voids 
occur next to condensations with baryonic clustering mass 



M ci --p B (d H ) - — 



1.7- 10 i4 Mr. 



0) 



(5) 



which corresponds to the baryonic mass of fat galaxy clus- 
ters (cl). The Reynolds number becomes 



dfiVs _ 97r 3 ^ 5 c 2 ag m d a H c V s n e p B 
v 5m 2 C(3) (k B Tf 



(6) 



At z v f it equals 158, somewhat above critical. At the 
boundaries of the clumps it is much smaller, R c (r, z) = 
R (z) [p B {f, z)/p B {z)] , where z codes the time, r the lo- 
cal position and the uniform terms refer to the would- 
be uniform state. While R = 158 at z v f is already not 
large, fragmentation leads to small values at the bound- 
aries, which enhances the effect. As the voids expand, 
baroclinic torques at their boundaries produce vorticity 
and turbulence due to misalignment of pressure gradients 
and density gradients. Pressure gradients will be normal 
to void boundaries, but density gradients need not. The 
rate of vorticity and turbulence production at the expand- 
ing protosupercluster boundaries is du/dt = Vp x Vp/p 2 . 
Observations of the Hubble Ultra Deep Field [17] show 
chains of protogalaxies and spiral clump clusters, as well 
as DM filaments, formed in this way. 

A connection to turbulence was established [18] in a 
study of the CMB temperature difference between two 
points at angular separation r, viz. (|AT| P ) ~ r^ p , where 
the average is taken over angles between 0.9° and 4°. For 
0.1 < p < 3 the exponent reads ( p w p/3, as in turbulence. 



A test of Gaussianity in CMB, (|AT|p) - (|AT| 3 )^ re- 
veals a marked deviation from the Gaussian value ( p = p/3 
in the interval 3 < p < 12, with £12 ~ 2.8, and coinciding 
with the C P of turbulence [19]. In rcf. [20] it is deduced 
that the data for the first CMB peak involve R ~ 100, in 
striking agreement with our estimate ~ 158. 

The pancaked structure of matter in between large voids 
arises dynamically since voids expand more than matter. 

The Axis of Evil. — WMAP data shown alignment 
of low order multipolcs (£ = 2,3) of the CMB spec- 
trum [3,21]; further coincidences occur up to t = 17. In 
our GHD approach the AE reflects density gradients of big 
bang turbulence [22] and mixing [23] subject to universal 
similarity laws of fossil turbulence and turbulent mixing 
[11]. Such turbulence fossils appear in CMB spherical har- 
monics axes and the spin axis of local galaxies clusters. 

Towards decoupling. — In the period near last scat- 
tering, Helium is already formed, so the density of protons 
plus H-atoms is n = 0.76 p B jmtq. The fractional ioniza- 
tion X = n e /n evolves according to Eq. (2.3.27) of [15], 



AX 
AT 



na_ X 2 -(l-X)/S f„_T An\ X_ 
HT 1 + A + { n AT T 



A = 



nA 3 e 157,894K/T ^ = % 



1.4377 10- 16 (T/K) 
1 + 5.085 lO-^T/tf) - 5300 

r 2s + 87riJ/[A 3 n(l-A)]' 



2tt 



m e k B T' 



(8) 



Here S is the Saha function and At the thermal length, 
while a and A are factors involving = 8.22458 s -1 the 
two-photon decay rate of the H 2s level and X a = 1215 A 
the Ly Q wavelength. We added the last term in (7) in order 
to allow that n ^ const. T 3 . Baryonic matter will expand 
less after clusters have formed. Let us take the geometric 
mean between no and full expansion, thus assuming that 
the matter clumps expand up till last scattering at zl by 
a factor yjaja^i < 2.2, implying p B — (a v f/a) 3 ^ 2 p B l . 

Initially S < 1, so the Saha law X = 1 - SX 2 con- 
tinues to hold. H-formation makes X decrease apprecia- 
bly, from where on we have to solve Eq. (7). The con- 
dition for maximal probability of last scattering [15] can 
be formulated as AJ/dT = J 2 , where J = carnX / HT in- 
volves the Thomson cross section ctt = Qir{Ti/ a em m e c) 2 = 
6.6525 10~ 29 m 2 . This fixes the surface of last scattering 
at T L = 2862 K, z L = 1050, compared to z dc = 1090 from 
WMAP5, and taking plac e at age t L = 408,000 yr. The 
clump size L c \ — d H c (z v f)y / aL/a v { corresponds to an angle 
c i = 180° L ci /ndA(zL) = 0.84°, or spherical index t c \ = 
180°/6» c i = 215, which agrees with the first CMB peak. At 
this moment X = 0.01 makes the Reynolds number as low 
as 0.12, thus exhibiting turbulence throughout the bulk, 
and predicting more CMB turbulence at smaller scales. 
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Magnitude of CMB temperature fluctuations. 

The smallness of CMB fluctuations, ST/T ~ 1CT 7 is one of 
the mysteries of cosmology. Indeed, how can it be consis- 
tent with a mass contrast of almost 100% between clumps 
and voids? Presently it is described by inflation, where 
its size is adjusted in the initial fluctuation spectrum [15]. 
In order to explain it from a physical mechanism, let us 
notice that not all clump energy can associate with tem- 
perature fluctuation, since in empty space the tempera- 
ture already decays with the rcdshift, T(z) = (z + 1)Tq. 
Compared to voids, extra energy of a clump that is avail- 
able for photons must stem from the kinetic energy of 
the protons, electrons and He atoms. Their densities are 
low, at decoupling ~ (1 + Zd c ) 3 n B p c /mN <~ 300/cc with 
p c = 1.04 • 10~ 26 kg/m 3 the critical density. At a temper- 
ature T < T v f this amounts to an excess kinetic energy 
of |fcs(T v f — T) for each of them, which corresponds to 
an energy density 2.37 £l B (p c /m N )k B T (l + z^ f )(z v{ - z). 
Let us assume that due to scatterings this gets redis- 
tributed to the local photons, thereby causing a pertur- 
bation in the photon energy density at last scattering , 
Su y {z L ) = An jPc c 2 (z L + l) 4 (ST/T) L . This yields 



ST 
T 



= 0.593 



n B k B T Q (zyt + l) 3 (z vf - z L ) 



Cl 7 m^c 2 



(ZL + I) 4 



(9) 



Though k B T n /m,]yc 2 = 2.5 • 10 13 is very small, we find 
for Zyf = 5120, z L = 1050 that 5T/T\ L = 6.1 • 10~ 8 , 
which corresponds presently to ST C \ = T c \ — T VO id = 
+0.17 /IK. This is the right order of magnitude, since 
5T 2 X = 204- [27r/4i(4i + 1)] pK 2 compares to the first peak 
of the correlator, C 22 o = 5800 • [2tt/4i(4i + 1)] pK 2 from 
WMAP5 [13, 16]. Voids do not have this baryonic con- 
tent, which explains the observed connection hot spots - 
(supcr)clusters, cold spots - voids. 

Fragmentation in the gas at two scales. At last 
scattering, the plasma turns into a neutral gas and further 



baryonic structures form. The free fall time is r„ 



1.68 



Myr, while the age is ft, = 0.41 Myr. The sound speed 
of a monoatomic gas is V s — y/5p/3p. For H with 24% 
weight of He, p = 0.82pk B T/rriN yields Vs — 5.68 km/s. 
The gas fragments at the Jeans scale Lj = V s r g = 9.78 pc 
into Proto-Globular Clusters (PGCs) with Jeans mass 



Mpgc — -PbL 
6 



6G 3 / 2 /o]/ 2 



38,000M Q . 



(10) 



This Jeans cluster formation is well known, but not always 
welcomed. In our approach it is a standard fragmentation. 

At decoupling the viscosity decreases from photon vis- 
cosity values to hot-gas values. The He viscosity can be 
estimated as r}a e (T L ) = 5.9 • 10~ 5 kg/ms. For the 76:24 
H-He mixture it will be about 0.76/8 + 0.24/4 = 0.155 of 
this. The critical viscous length Lsv = (VsV/Gp 2 B ) 1 ^ 3 = 
3.9 • 10 14 m implies a condensation mass 



M = ^ PB L SV = = 13M ffi = 3.9 • 1O" 5 M . (11) 

We may call these objects H-Hc planets, Primordial Fog 
Particles or Milli Brown Dwarfs. Their mass is in good 
agreement with estimates from microlensing of a distant 
quasar [24, 25] and so-called comctary knots in the Helix 
nebula [26-28]. It was anticipated both by Gibson [10] 
and Schild [24] that galactic dark matter is composed from 
such planets. Each PGC contains about a billion of them. 

Galaxies. — We may relate galaxies to the Jeans 
mechanism at the end of the plasma epoch. The sole rel- 
evant aspect is then the decrease of the speed of sound 
from plasma values to hot gas values. Taking the geomet- 
ric mean velocity V s = {Vg lasma V§ as ) 1 / 2 = 874 km / Sj we 
get the Jeans scale Lj = 1.5 kpc and corresponding mass 



M, 



gal 



6 PbL > 



6G 3 / 2 Ps /2 



1.4- lO n M . 



(12) 



x We neglect here that H atoms have formed in mean time. 



The corresponding CMB angle is 0q = 4.7' and the angu- 
lar index is Iq = 2300. For this mass regime the formation 
time is limited, because the sound speed continues to de- 
crease to gas values, from where on PGCs are formed. This 
explains why a lot of baryons are not locked up in galaxies 
with their baryonic dark matter, but located in intraclus- 
ter and intercluster X-ray gas. That gas has become hot, 
with temperatures up to 100 keV (1 keV/k B = 1.16T0 7 K), 
due to virialization after neutrino condensation on clusters 
at z ~ 28 or t vc = 120 Myr [8]. At such high temperatures 
the gas may allow nuclear fusion up to tellurium [29]. 

Role of PGCs. — In some of them, still warm, colli- 
sion processes have quickly led to star formation, basically 
without a dark period, thus transforming them into OGCs. 
Other PGCs transformed in ordinary stars. In the major 
part of the PGCs the planets have frozen and still per- 
sist without stars. These PGCs are in virial equilibrium 
and act as ideal gas particles that constitute the galactic 
dark matter. Their physical presence explains why the 
isothermal model describes the basic features of galactic 
rotation curves so well, that is, linear growth at small ra- 
dius, plateau at large radii. To improve the fit, one may 
consider mixtures with several isothermal components [8] . 

In the centers of galaxies near passings of PCGs will 
heat their planets and induce star formation. Since this 
is mainly a two-particle effect, the luminosity of a galaxy 
is expected to relate to the PGC mass density as I ~ 
Jd 3 

r Ppcg- m the isothermal model ppec — c 2 /27rGr 2 , 
where a v is the velocity dispersion, so the 1/r 4 fall off of 
the integrand makes the luminosity finite. This results in 
the scaling L ~ tr 4 /i?, which is the Faber- Jackson relation 
[30], with an additional characteristic bulge scale R. The 
Tully-Fisher relation is likewise explained, as it involves 
the rotation velocity, which scales with a v . 
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In several instances the matrix of dark PGCs is revealed 
by new star formation. When agitated by tidal forces 
the collision frequency of the planets will increase caus- 
ing re-evaporation of the frozen gases, increased size and 
friction, and the possibility of planet mergers to produce 
larger planets and eventually new stars. The existence of 
galaxy dark matter in the form of clumps of frozen primor- 
dial planets is clearly revealed in photographs of galaxy 
mergers such as Tadpole, Mice and Antennae. On the 
respective photographs one can see numerous bright clus- 
ters of comparable size, that are identified here as PGCs 
turned into young globular clusters. They are located in 
star wakes as the merging galaxies enter each others dark 
matter halos and heat up the planets in the PGCs on 
their path through the dark matrix. The effect exists only 
within a certain radius, the boundary of the PGC cloud. 

Role of H-He planets. — From the GHD scenario 
following decoupling, the first stars form gently by a fric- 
tional binary accretion of still warm PFPs to form larger 
planet pairs and finally small stars as observed in OGCs. 
Thereby they create an Oort cavity as clearly exposed in 
e. g. the Helix planetary nebula. Slow turbulent mixing 
from the rain of planets will not mix away the dense car- 
bon core. By conservation of angular momentum the star 
spins rapidly as it compresses producing strong spin radi- 
ation along the spin axis and at the star equator. Stars 
can form as binaries in PGC clumps. 

At decoupling the entire baryonic universe turns to a 
fog of H-He planet mass clouds. Due to the expansion 
they cool and the freezing temperature of hydrogen and 
helium occurs at redshift z w 30 a few hundred million 
years later, producing the frozen dark baryonic planets in 
clumps predicted as the galaxy dark matter [10]. Neu- 
trino condensation at z ~ 28 coincides with this, showing 
that the extra-galactic, gaseous H-He planets heat up to 
become hot X-ray gas. That the dark matter of galaxies 
should be planets of earth mass was independently pro- 
posed by the Schild 1996 interpretation of 5 hr twinkling 
periods in quasar microlensing observations [24,25]. Thou- 
sands of these planet crossings have been observed by now. 

As the universe expands the planets that did not turn 
into stars freeze and the PGCs become less collisional and 
diffuse out of the 3 kpc scale protogalaxies to form the ob- 
served typical R = 100 kpc dark matter halos with a basi- 
cally isothermal distribution. Consistency of this picture 
is shown by the isothermal estimate M ga \ = 2cr%R/G ~ 
1.9 • 10 12 M Q for o v = 200 km/s. 

Spiral galaxies reflect the accretion disks of dark halo 
PGCs frictionally spiraling back toward the Lm scale pro- 
togalaxy remnant of ancient stars [31]. Violent massive 
stars form at Lst = e 1/>2 / (pG) 3 ^ 4 as turbulent mael- 
stroms at protogalaxy cores giving spin radiation quasar 
and gamma ray burst events. For e <~ 10 -10 m 2 /s 3 one 
derives a billion solar mass quasar core. 

When frozen, the H-Hc planets are too small to dim 
light, even from remote sources, but they can account for 



dimming when they are heated. Warm atmosphere diame- 
ters are w 10 13 m, the size of the solar system out to Pluto, 
bringing them out of the dark. The separation distance 
between planets is w 10 14 m , as expected if the PGC 
density of planets is the primordial density p = 2 x 10~ 17 
kg m~ 3 . In planetary nebula such as such as the nearby 
Helix, dark planets at the boundary of the 3 x 10 15 m 
Oort cavity are evaporated. HST optical images of Helix 
show w 10 4 "cometary knots" , planet-atmospheres w 10 13 
m which we identify as H-He planets with metallicity, and 
Spitzer shows > 10 5 in the infrared from the 10 3 Mq avail- 
able [12]. Meaburn et al. determine for cometary knots in 
Helix a mass of 1. 1O _5 M , and conclude that the glob- 
ules and tails are dusty [26], and later report a mass of 
2. 10 _5 M Q [27]. Huggins et al [28] measure a mass of 
5. 10~ 6 M Q via the radio measurement of CO emission. 
These findings support the GHD prediction of earth mass 
planets as the baryonic dark matter repository. 

Dimming by dense 10 13 m planet atmosphere gases or 
realistic quantities of dust is negligible at the 20% levels 
observed [32]. Such large dimming of obscured lines of 
sight observed in a planetary nebula (e.g. Helix) requires 
post turbulent electron density forward scattering [33] as 
observed by radio telescope pulsar scintillation spectra, 
embedded in the Kolmogorovian "great power law on the 
sky" [11], which can now be understood from GHD as 
remnants of supernova powered turbulent mixing in our 
local PGC [33]. The planet atmosphere cross section for 
SNe la dimming a w 10 26 m 2 gives a photon mean free 
path 1 j tig ~ 3 x 10 15 m from the primordial PFP number 
density n , comparable to the observed Helix planetary 
nebula shell thickness and consistent with w 5% of the SNe 
la lines of sight unobscured. Clouds of warm H-He planets 
can be responsible for the Ly Q forest - not hydrogen clouds 
that have remained undetected. 

Conclusion. — Within gravitational hydrodynamics 
(GHD) for a flat space Friedman cosmology with (hot) 
neutrino cluster dark matter we have described various 
structures that form due to hydrodynamic instabilities in 
the baryonic plasma. The approach accounts for a wealth 
of observations and relations between them. 

The first step, when photon scattering from free elec- 
trons makes the plasma too viscous, is the creation of 
proto-voids at z v f = 5120, before decoupling, that ex- 
pand to present average cosmic voids. Proto-supervoids 
occur too, but are rare. The matter condenses in proto- 
galaxy-clusters while large scale turbulence at boundaries 
of expanding voids in the plasma respect the earlier de- 
termined axis of evil. The coincidence of turbulence 
properties in the CMB and turbulent fluids supports the 
GHD model. Pancaked galaxy (cluster) structures in be- 
tween large voids arise dynamically, since the voids expand 
faster, while the matter remains gravitationally bound. 

The assumption that the kinetic energy of the protons, 
He atoms and electrons is redistributed to the photons, ex- 
plains the micro-Kelvin scale of CMB fluctuations and the 
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connection cold spots - voids, hot spots - (super)clusters. 

Before the viscous instability the plasma was quite 
homogeneous, due to the free streaming neutrinos that 
damped inhomogeneities at the free streaming scale of 10 
kpc at z v f. So, irrespective of their location on the sky, 
the fluctuations caused by the viscous instability have the 
same order of magnitude. This well known large angle 
correlation is in GHD an effect of simultaneity (of void 
formation in the homogeneous plasma), not of causality. 

At decoupling first galaxies are formed and later primor- 
dial globular clusters (PGCs) of about 38,000 solar masses. 
Some of them turn into old globular clusters (OGCs) and 
others form the stars in galaxy bulges. But most PGCs, 
millions per galaxy, remain dark and constitute the galac- 
tic dark matter as an ideal gas, which explains why the 
isothermal model describes galactic rotation curves well. 
A galactic dark PGC matrix also accounts for numerous 
young globular clusters seen in galaxy mergers. Near PGC 
crossings in the center of galaxies will form stars, which 
explains the Tully-Fishcr and Faber- Jackson relations. 

At the decoupling the viscous scale becomes much 
smaller which turns all matter into H-He planets of a few 
earth masses. Some of them coalesce to form the first small 
stars, but most freeze to earth scale. The MACHO [34] 
and EROS [35] collaborations have searched in vain for 
such objects. Still, they are not excluded because they 
do not occur uniformly but in PGC clumps. Theoretical 
descriptions of clumped MACHOs in the dark halo were 
started in [36]. But thousands of planets have been ob- 
served in microlensing [24, 25] and, reheated, thousands 
more in planetary nebula such as Helix [26-28,37]. Hot 
H-Hc planet atmospheres may dim distant supernovas. 

We have adopted one set of cosmological parameters, 
which performed rather well, but not attempted an op- 
timization. While in the CDM model the main cause of 
clustering is dark matter with baryons a second order ef- 
fect, the GHD scales will be rather sensitive to the precise 
parameter values. Large scale numerical hydrodynamics 
simulations of separate steps of the fragmentation process 
are expected to result in precise fits for the mass fractions 
of baryons and neutrinos, and the Hubble constant. 

We have not considered the large scale power spectrum 
[15]. We may recall that there are reasons to question 
whether baryons trace the neutrino dark matter well [8]. 

Let us finally see how some problems of the 
ACDM paradigm mentioned in the introduction are solved 
naturally in GHD. Population III stars may have been 
rare, since reionization may find its origin in neutrino con- 
densation on galaxy clusters [8] . Dwarf galaxies with a lot 
of (baryonic!) dark matter may pertain to PGCs with 
incomplete star formation. The related fact that OGCs 
often exhibit stars formed at several epochs is likewise 
explained by further sets of reheated, pre-existing H-Hc 
planets. Correlations in galaxy structures are expected 
since they all form early; baryon acoustic oscillations do 
not show up in GHD. Galaxy formation when the universe 
was 4-5 billion years young may refer to proto-galaxies 



with late-stage star formation by close PGC encounters. 
Dwarf satellites that swarm our galaxy just like its stars 
may just relate to single PGCs with modest star forma- 
tion, popping out of the matrix of dark PGCs. 
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